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Abstract 

The information content of a spinning sound field is analyzed using a combination of exact and 
asymptotic results, in order to set limits on how accurately source identification can be carried out. 
Using a transformation of the circular source to an exactly equivalent set of line source modes, given 
by Chebyshev polynomials, it is found that the line source modes of order greater than the source 
wavenumber generate exponentially small fields. Asymptotic analysis shows that the remaining, lower 
order, modes radiate efficiently only into a region around the source plane, with this region shrinking 
as the mode order is increased. The results explain the ill-conditioning of source identification methods; 
the successful use of low order models in active noise control; and the low radiation efficiency of subsonic 
jets. 

1 Introduction 

Source identification, the problem of determining an acoustic source from field measurements, has been 
attempted using a number of approaches in various different technologies. This paper examines the problem 
of identifying the source which generates a spinning acoustic field. The source which generates such fields 
can be represented as a set of modes which vary with azimuth on a circular disk, whether or not the 
source includes a spinning element. Examples include rotating systems such as cooling fans [Tl[2], helicopter 
rotors fSl*?] duct terminations such as aircraft engine intakes O [6l [71 [H [9l [10] and jets [TT] if a jet is modelled 
distribution of disk-shaped sources. 

There are two broad categories of problem where source identification is required, corresponding to 
'forward' and 'backward' projection of the field. In the forward problem, the aim is to estimate the source 
distribution accurately enough to allow the field to be predicted at positions other than the original measure- 
ment points. This has been done, for example, in extracting source parameters from near-field measurements 
of propeller noise, with the parameters then being used to calculate far-field noise[12'. 

In backward projection the problem is determination of the source proper from field measurements. This 
might be done in order to decide on noise control measures or because the acoustic source corresponds to 
some other physical variable of interest. In the first case, the aim is usually to find the acoustic source 
strength distributed over the source region, such as a rotor diskt]jL2j or the termination of a duct[71[51[51ll0j. 
An example of the second application is the study of noise generation by turbulent iets|ll). where the aim 
is to determine the fluid-dynamical mechanisms which give rise to the acoustic source. 

Iir any case, it is well knowir that the problem of source identification is poorly-conditioned, meaning 
that small measurement errors can give rise to very large changes in the estimated source. Previous studies 
of the structure of spinning fields[13[ [T^ [T51 [T5] have shown that the field decays exponentially away from 
the source. This means that in the forward projection problem, errors in the estimated acoustic source will 
decay and the predicted acoustic field may well be quite accurate, even if the source is not well recovered. 
On the other hand, in the backward projection problem, the exponential decay moving away from the source 
corresponds to exponential growth moving towards it, leading to large errors in the estimated source. 



*This paper is based on work which was to have been presented at the 159th Meeting of the Acoustical Society of America, 
Baltimore, 2010. 
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Figure 1: Coordinate system for radiation prediction 



In this paper, the radiated field from a disk source is analyzed to examine how much information 
about the source can actually be detected in the field. In a previous paper on a possible method for source 
identificationflT] . it was shown that the far- field noise is band-limited Fourier transform of a line source 
which is exactly equivalent to the disk source. In this paper, without recourse to far field approximations, it 
is possible to establish fundamental limits on the number of degrees of freedom of a field, limits which are 
determined by the source frequency. The implications of the analysis are discussed with respect to some real 
source identification and radiation prediction problems. 



2 Analysis 

The problem considered is that of the field generated by an azimuthally varying distribution of monopoles 
with strength s{ri,di) given by the Rayleigh integral [TBI ^13] 

p{r,e,z,uj)^ I I s(ri,0i)— d^indn, (1) 



JO Jo 

= r^ +rl- 2rri cos(6' - 6*1) + z^, 

where the source is distributed over the unit disk in the plane z = 0, variables of integration have subscript 
1 and the coordinate system is shown in Figure [T] The wavenumber k = uj/c, and c is the speed of sound. 

Taking one azimuthal mode of the source distribution, s(ri,0i) = s„(ri) exp jn^i, the radiated field for 
one mode can be written p = p„ exp jnfl: 

Pn(?',2) = y J s„(ri) — ^-^ — d6'iridri, (2) 

= r'^ +rl- 2rri cos 6*1+2^. 

The integral of Equation [3] has been extensively studied due to its relevance to rotor acoustics and, un- 
der suitable conditions, as a good approximation to radiation from ducts. Many problems in source 
identification[TJ [21 [H [H [7l [8j [9j [10] can be viewed as attempts to recover the source term s„(ri) from 
measurements of p„ . 

The remainder of this paper consists of an analysis of the integral of Equation[2l which will establish limits 
on the information about the source which is radiated into the acoustic field, thereby fixing how accurately 
a source can be identified. The results are also applicable to the question of the detail with which a source 
need be specified in order to accurately predict the acoustic field, and to that of the radiation efficiency of 
jets. 
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Figure 2: Coordinates for transformation to equivalent line source 



2.1 Equivalent line source 

The first step in the analysis is to transform the disk source into a line source which generates (exactly) the 
same acoustic field. This is a transformation which has been used, in the axisymmetric case, in studies of 
transient radiation from pistons[19' '20], and, with azimuthal variation, in studies of rotor acoustics[131 1141 
[TS1[T2]. The first stage is to switch from the source-centred cylindrical coordinates {r,9,z) of Figure [T] to 
the observer-centred coordinates {r2,02,z) of Figure[21 Under this transformation. Equation [2] becomes, for 
r > 1: 

Pnik,r,z)^ I '—^K{r,r2)r2dr2, (3) 



f-r+l 




Jr-1 








{rl + 




-1 




47r Je 


e 

(0) 



'^Sn(ri)d^2, (4) 

where the source function K(r,r2) depends on r, the observer lateral separation, but is independent of z, 
the axial displacement. The coordinate systems are related by: 

rj ^ r'^ + +2rr2Cos92, (5a) 
1 To sin 02 

01 = tan^i V- 5b 

r + r2 cos 02 

and the limits of integration in Equation 2] are given by setting ri — 1: 

0T=co.-^'-^. (6) 
The function K[r, r2) has square-root behavior at its end-points [17J, r2 — r ± 1, and can be expanded: 

oo 

K{r,r2) = u,„(r){/„(s)(l - s^f'^ (7) 
with s = r2 — r and Um{s) the Chebyshev polynomial of the second kind. 
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Inserting this expansion into Equation [3] 



p„(fc,r, z) = ^ u,n{r)Irn{k,r,z), (8) 

m=0 

I^{k,r,z)^ j ^^—Ura{s){r + s){l~s'Y'^ds, (9) 

Ti- gjfc_R 

— — (r + cos;9) sin(TO + 1)/? sin/3d,3, (10) 
it 

= (r + cos/3)2 + 

with the transformation s = cos /? and use of the definition of the Chebyshev polvnomial|21). Um{s) = 
sin[(m + sin/3. 



3 Radiated field 

The analysis of the previous section gives us a model of a spinning acoustic field expressed in terms of 
an exactly equivalent line source composed of a superposition of modes given as Chebyshev polynomials, 
with the modal coefficients functions of observer radius r. but not of axial displacement z. In this section, 
we use the model to draw basic conclusions about the acoustic information which is available for source 
identification. 



3.1 Cut-off modes 

The first conclusion we can draw from the integral expression for is that modes with m > k are exponen- 
tially small and can be considered 'cut-off'. For z = 0, Im{k, r, 0) can be evaluated exactly using EguationfTUl 
and standard relations for Bessel functions|21): 

/,„(fc,r,0)=r^(™-f l):^^^^ej^'-. (11) 

For m + 1 large and k < m + 1, the Bessel function Jm+i{k) decays exponentially, i.e. the higher modes 
are 'cut off'. Since llmif-, z)\ has its maximum in the plane z = 0, we can further conclude that modes with 
m > k are cut off everywhere and cannot be detected in the field. This is an exact result which places a first 
limit on the information radiated. 



3.2 Cut-on modes 

A second limit on the information available in the acoustic field can be found by asymptotic analysis of Im 
which can be rewritten: 

Irn = (Qm+2 (k, r, z) + Q ^,n-2ik, T, z) 

- Qm{k,r,z) - (3_„(fc,r, z))/4, 

with: 

Q™(fc,r,z)=^^d'=^(«!:±^d/?, (12) 

m = i? + 7/3, 
7 = m/k. 
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Figure 3: Cut-ofF lines for 7 = to/8, m — 1, 2, . . . , 6, 7; exact solution solid; Zc = ar/j dashed. 



The integral Qm is in a suitable form for stationary phase analysis [551, which depends on finding the station- 
ary points of ip, i.e. values of /? where dip/dP = with < /3 < tt. Upon differentiation and rearrangement, 
the condition dijj/df3 — takes the form of a quartic equation: 



(13) 



where C — cos/3 and = 1 — j'^. To lie in the domain of integration, the stationary phase points must be 
real with |C| < 1. This leads to the requirement that < 7 < 1 and \z\ < Zc, a 'cut-off' value of observer 
axial displacement beyond which the phase function ip has no valid stationary points. The two values of C, 
denoted C+ and C_, at the limits oi \z\ are: 



C±(z) 



zba, 

-r/4+ (r2 + 8^2)1/2/4, 



z 



-0, 

= Zc- 



Denoting Cc 



C+(zc) — C_(zc), the cut-off value of z is: 



Zc = (a- -U~r'~(r + Cc)/j, 
ar/"f, a/r 0. 



(14) 



(15) 



sin ^ 7. 



Written in spherical coordinates, the asymptotic cut-off lines Zc — ar/j are rays with polar angle ' 
For completeness, we note that for 7 = 6, there is no cutoff and the line source mode radiates into the whole 
field with amplitude proportional to A:"^/^ 

Figure [3] shows the cut-off lines, exact and asymptotic, for < 7 < 1. The accuracy of the asymptotic 
approximation for the cut-off line is confirmed and the plot shows which radiated modes are detectable in a 
given part of the acoustic field. 

Using the previous results, the asymptotic behavior of the basic integral is given by: 



•3/2_ 



+ j 



1/2. 



{kRlyT^ 



27r 



27r 



1/2 



(r -I- C+) 



C-{r + 2C-) 
, ip±^ R±- 



1/2 



7Cos~^ C±, 



(r + C_), fc -> 00, 



(16) 



and: 



{Qm+2{k,r,z) 



,{k,r,z))/A, 



(17) 



where Q-m and Q-m-2 are neglected since they have no stationary phase points and decay much faster 
than Qrn and Qm+2 with increasing k. 
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Figure 4: Numerical and asymptotic evaluation of /„, r = 1.125, fc = 8, n = 1: left hand side real, right 
hand side imaginary; solid line numerical evaluation, dashed line asymptotic. Equation 1171 



Figure O compares the asymptotic approximation for /i(fc, r, z) to a numerical evaluation of the integral. 
The cut-off value of z for 7 = 3/fc is indicated and the real and imaginary parts of the integrals are plotted 
separately. As z — > Zc, the stationary phase approximation to Qm+2 breaks down and there is a resulting 
loss of accuracy. Away from this point, however, the approximation to is accurate on both sides of Zc- 

The asymptotic analysis shows that the cut-on line source modes, those with m < k, radiate efhciently 
into a region bounded by ±Zc(r, 7). This gives a second limit on the information available in the acoustic 
field. 



3.3 Far field approximation 



For completeness, we give a far field approximation of the line source radiation integral, valid outside the 
region covered by the asymptotic expansions of ^3.21 Using the standard approximation, R ^ Rq + [r 
r2)sin0, l/i?« 1/Ro with i?g 



j,2 p. 



1^1 



1 



i J TT — : — 

Rq k sin ( 

- jJm(A:sin0) 



r+i 



im + 2) 
ksind) 



Jm+i(fcsin 4>) 



(18) 



so that in the far field, Im decays as k ^. 



4 Information in spinning sound fields 

Summarizing the results of the previous section, the nature of a spinning sound field is seen to be determined 
by its wavenumber k and the relation of k to the set of modes contained in the equivalent line source. The first 
result, that line modes with m > k generate exponentially small fields, means that the acoustic field contains 
a limited amount of information about the source. Such a result has been derived previously by showing 
that the far field pressure is a band- limited Fourier transform of the line source strength |17|. but this new 
result establishes an exact limit on the information in the field, without needing a far-field approximation. 

The second result, from the stationary phase analysis, shows that the modes which do radiate, those with 
m < fc, are more efficient in some parts of the field than in others. The higher order radiating modes are 
detectable only near the source plane, with a lower radiation efficiency at larger z. The only mode which 
radiates efficiently over the whole field is the 'plane' mode m = 0. 

The following sections discuss some implications of these findings for different problems. 
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Figure 5: Line source reconstructed with different modes for r — 5/4: tliick line full source term K{r, r + s); 
thin lines K{r,r + s) reconstructed using 1-5 modes. 

4.1 Low speed rotors 

One result which is of immediate interest is that, in some sense, low speed rotors have the same acoustic field. 
Given that for a system of radius a rotating at angular velocity J7 the non-dimensional wavenumber of the 
nth harmonic of the radiated field k = nVLa/c — nMt, with Mt the source tip Mach number, low speed rotors 
will have fc < 1 over the first few harmonics of the signal. This means that the sound field is dominated by the 
zero order line mode and any set of rotors of a given blade number operating at the same speed, whatever 
their blade geometry, will have the same acoustic field, to within a scaling factor. Indeed, experiments 
in active noise control of noise from cooling fans [23] have found good results by discretizing the inverse 
model of the fan into three sections, i.e. using three degrees of freedom, for a value of ka ~ 0.8. The same 
group, in an earlier study of the conditioning of the inversion problem, using a hemispherical arrangement of 
microphones [T] found that the condition number reduced as frequency was increased, a finding they ascribe 
to "insufficient spatial resolution of the source for frequency below 200Hz." In the light of the analysis above, 
an alternative viewpoint might be that as frequency is increased, more line source modes are cut on and the 
information available in the acoustic field becomes more nearly sufficient for source reconstruction. 

4.2 Source identification 

The original motivation for this work was the problem of identifying a rotating source. The results of ^and 
^ can be used to help show why this is a hard problem and to indicate how it might best be approached. 

The first obvious consequence of the result of ^3.11 is that the acoustic field has a limited number of 
degrees of freedom. For a field of given wavenumber k, no more than k modes can be detected in the field, 
i.e. the field has M degrees of freedom, with M the largest integer M < k. Attempting to identify sources 
using more than M degrees of freedom is inherently ill-conditioned because the modes with m > k generate 
exponentially small fields. 

Secondly, the asymptotic analysis shows that much of the information in the field is not detected by 
microphones. For a given microphone polar angle, only a subset of the M radiating modes will be easily 
detected. If the microphone is at polar angle </>, modes with 7 > sin (f> may not be detected or, alternatively, 
the microphone only detects radiation from modes of order m < k/ sm(f). 

The effect this has on the source identification problem is shown in Figure [5] This shows, as a thick 
line, K{r,r + s) for r = 5/4, n = 2 and s(ri) = 1, calculated using previously published closed-form 
expressions [TJ. The thin lines show K reconstructed using the first M terms of the sum in Equation [71 with 
M = 0, 1,2, 3, 4. The convergence towards the exact value of K is obvious, but it is also obvious that this 
convergence is so slow that a large number of terms will be needed in order to accurately reconstruct K. 
Given that the number of modes which can be detected depends on wavenumber fc, it is clear that except at 
very high frequencies, only low resolution source reconstruction is possible. 
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Figure 6: Acoustic integral with exact line source (solid line) and lower order modes (dashed lines) for 
m < 0, 2, 4, fc = 2.5, r — 1.25; real part on left, imaginary part on right. 
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Figure 7: Acoustic integral at 2 against maximum mode order M included; real part solid, imaginary 
part dashed. 




4.3 Source resolution 

The observations of the previous subsection give a pessimistic outlook for source identification: the informa- 
tion necessary for accurate reconstruction of a source is simply not available, even in theory. On the other 
hand, if the objective is to reproduce an acoustic field, for noise control, say, only limited knowledge of the 
source is necessary. Looking again at Figure [SJ although there is a large difference between the exact line 
source and the source produced by summation of lower order modes, the acoustic field generated by the 
lower order modes with m < k will be indistinguishable from that generated by the exact source, since the 
modes with m > fc do not contribute to the radiated sound. 

Figure [6] shows the acoustic field on a sideline r = 1.25 computed using the exact line source of Figure [5] 
and by summation of the field due to the lower order modes, with m < 0,2,4 shown. The wavenumber 
k = 2.5, so that cutoff begins at m = 3. Indeed, the sum over the first three modes is very close to the exact 
result, and when the first five modes are included, the result is indistinguishable from the exact field. 

Figure [7] shows the development of p„ at z = as more modes are included in the summation. It is 
clear that the sum has all but converged after the m = 2 mode is added and is practically unchanged after 
the m = 4 mode is added, confirming the conclusion that the higher order modes are cut-off and do not 
contribute to the acoustic field, in spite of their quite high amplitudes. The clear conclusion is that the 
acoustic field depends on the lower order modes and sources which differ only in the higher order terms of 
their line source decomposition cannot be distinguished by a source identification procedure. The corollary 
of this statement is that for an accurate prediction of the field, sources need only be specified to a resolution 
sufficient to correctly identify the amplitudes of the lower order modes. 
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4.4 Jet noise 

In a recent paper [TT]. Jordan et al. examine noise production by a turbulent jet using proper orthogonal 
decomposition (POD) to perform a modal decomposition of the flow and a technique called MOD ('most 
observable decomposition') to decompose the acoustic far field. They find that while 350 flow modes are 
needed to account for 50% of the turbulent kinetic energy in the flow, only 24 modes are needed to account 
for 90% of the acoustic energy. If the jet is viewed as a distribution of disk sources along the jet axis, 
the results of this paper show that we should expect that only a small fraction of the modes will radiate 
noise and that in a complex source such as a jet, the bulk of the modes will have m > k and will generate 
exponentially small fields. In a study of noise sources in a jet, Freundplj filters the source terms to leave 
"a set of modes capable of radiating to the far field", based on a wavenumber criterion, but he notes that 
"additional cancellation may occur due to the radial structure of the source". The analysis of the previous 
sections offers an approach for the identification of the radial terms which will give such cancellation. 

5 Conclusions 

An analysis of the information content of a spinning sound field has been presented. It has been shown, 
on the basis of an exact analysis, that the acoustic field around a spinning source has at most M degrees 
of freedom with AI < k, the acoustic wavenumber. This result arises from the replacement of the disk 
source with an exactly equivalent line source given by a sum of modes. Most of these modes generate 
exponentially small acoustic fields, i.e. are cut off, and the remaining, lower order, modes radiate efficiently 
only into sectors of the acoustic field which become smaller as the modal order increases. The results explain 
a number of features which have been observed in the literature, including: the possibility of using low order 
source models for noise control; the ill-conditioning of source identification methods; and the low radiation 
efficiency of subsonic jets. 
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